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I.  INTRODUCTION 

Coronal  magnetic  fields  are  intimately  linked  with  most 
aspects  of  solar  activity.  Open  field  configurations  are 
causally  related  to  coronal  holes;  strong  fields  are 
associated  with  solar  active  regions;  unstable  coronal 
fields  result  in  solar  flares.  Observations  of  coronal 
fields,  however,  have  been  based  largely  on  soft  X-ray 
images  such  as  were  taken  on  Skylab,  which  clearly  showed 
the  important  role  of  magnetic  loops  in  the  physics  of  the 
solar  corona.  However,  such  images  cannot  measure  the 
strength  or  the  fields.  Knowledge  of  the  field  strength  and 
distribution  is  important  since  in  solar  active  regions,  the 
magnetic  energy  density  greatly  exceeds  the  plasma  ("nkT") 
energy  density. 

Although  indirect  estimates  of  the  coronal  field 
strengths  have  been  oDtained  by  extrapolating  from 
photospheric  magnetographs,  uncertainties  in  calibration  and 
well-founded  misgivings  about  the  choice  of  appropriate 
magnetic  field  model  limit  the  usefulness  of  this  strategy. 
Microwave  observations  provide  another  approach.  Under 
normal  circumstances,  the  solar  corona  is  optically  thin  at 
microwave  frequencies,  resulting  in  observed  brightness 
temperatures  or  10*  to  10s  K.  In  the  presence  of  strong 
magnetic  fields,  however,  gyroresonance  opacity  can  render 
the  corona  optically  thick  so  that  coronal  brightness 
temperatures  are  achieved  at  such  locations.  The 
corresponding  sources  have  been  ooserved  by  many  workers 
with  high  resolution  microwave  imaging  instruments  such  as 
the  Very  Large  Array  and  the  Westerbork  interferometer. 

Gyroresonance  opacity  is  significant  at  frequencies 
which  are  low  integral  multiples  of  the  local  gyrof requency. 
It  is  important  to  recall  that  the  gyrof requency , 

[GHz]  «  2.8  B  [kg]  (1) 

depends  only  on  the  magnetic  field,  B,  and  is  independent  of 
electron  temperature,  density  or  velocity.  At  temperatures 
and  densities  typical  of  coronal  conditions,  gyroresonance 


opacity  usually  provides  significant  optical  depth  at  the 
2nd  and  3rd  harmonics  in  the  extraordinary  propagation  mode 
(left  circular  polarization  for  downward  directed  fields) 
and  at  the  2nd  harmonic  in  the  ordinary  mode.  The  optical 
depth  at  higher  harmonics  is  usually  much  less  than  1  and  so 
will  be  neglected  here. 

These  considerations  lead  to  the  schematic  picture 
illustrated  in  Figure  1,  whereby  at  a  fixed  frequency,  the 
coronal  microwave  emission  above  a  sunspot  is  generated  in 
magnetic  shells,  representing  isogauss  surfaces.  At  5  GHz 
for  example,  the  emission  would  come  from  shells  where  the 
magnetic  fields  are  approximately  600  and  900  gauss.  (At 
600  and  900  gauss,  the  gyrof requencies  are  1.7  and  2.5  GHz, 
whose  2nd  and  3rd  harmonics  are  5  GHz.)  Thus  at  5  GHz,  the 
ordinary  mode  emission  comes  from  a  surface  where  the 
magnetic  fields  is  900  gauss,  whereas  the  extraordinary  mode 
emission  comes  from  a  600  gauss  shell.  (In  the 
extraordinary  mode,  the  900  gauss  shell  also  would  emit  at 
5  GHz,  but  such  emission  normally  would  be  hidden  under  the 
opaque,  600-gauss  shell.) 

To  imaging  observations,  these  shells  appear  as 
localized  areas  coronal  temperatures.  Their  diameters  are 
determined  by  the  intersection  of  the  isogauss  shells  with 
the  transition  zone.  (Emission  by  the  shell  below  the 
transition  zone  is  unimportant  since  even  if  it  were 
optically  thick,  it  could  not  have  10*  K  brightness 
temperature.)  Note  that  at  higher  frequencies,  the  shells 
would  correspond  to  stronger  fields  and  so  we  would  expect  a 
decreased  diameter.  At  the  "transition"  frequency 
corresponding  to  the  snell  which  is  tangent  to  the 
transition  zone,  we  would  expect  that  the  diameter  would 
approach  zero  and  at  frequencies  above  this  cutoff,  coronal 
brightness  temperatures  would  not  be  achieved. 

The  variation  or  the  diameter  or  these  shells  with 
frequency  depends  on  the  configuration  of  the  magnetic  field 
in  the  corona.  Thus  measurement  of  the  diameter  and 
temperature  or  these  shells  as  a  function  of  frequency  can 
provide  a  direct,  quantitative  measurement  of  the  magnetic 
field  distribution  in  the  corona.  More  detailed 
calculations  illustrating  these  concepts  are  presented  in 
such  papers  as  Kruger,  Hildebrandt  and  Furstenberg  (1985). 

This  report  outlines  the  initial  results  of 
measurements  designed  to  measure  the  size  and  brightness 
temperature  or  these  shells  as  a  function  of  frequency. 


II.  OBSERVATIONS 


The  observations  were  made  with  the  frequency-agile 
interferometer  or  the  Owens  Valley  Radio  Observatory  (OVRO) . 
This  system  (Hurford,  Read  and  zirin  1984)  can  observe  at  up 
to  86  frequencies  between  1  and  18  GHz  in  both  right  and 
left  circular  polarization.  By  rapidly  switching  among 
frequencies,  detailed  spectra  can  be  obtained  in  a  few 
seconds. 

Early  results  (Hurford,  Gary  and  Garrett  1985)  used  two 
antennas  for  measurements  of  active  regions.  With  only  a 
single  baseline,  it  was  not  possible  to  directly  measure  the 
size  or  the  sources.  The  resulting  spectra  were  consistent, 
however,  with  the  arguments  advanced  in  Section  I.  Soviet 
workers  (cf  Akhmedov  et  ai.  1983)  have  conducted 
corresponding  multifrequency  measurements  between  7.5  and  15 
GHz  using  their  RAT AN-6 00  instrument. 

The  present  results  are  based  on  the  first 
interferometric  measurements  using  three  antennas.  As 
discussed  below,  the  resulting  three  baselines  permitting  a 
direct  measurement  of  source  size.  The  east-west  separation 
of  the  antennas  was  0.183,  1.067  and  1.250  km  with  the 
central  antenna  displaced  0.032  km  north  of  the  line  joining 
the  outer  two.  The  resulting  resolution  (fringe  spacing) 
was  inversely  proportional  to  frequency  (with  a  value  of  10 
arcsec  at  5  Ghz) ,  and  was  sufficient  to  resolve  the  sources 
at  each  frequency.  Although  data  were  acquired  at  86 
frequencies  between  1  and  18  GHz,  but  the  analysis  discussed 
here  will  be  limited  to  56  frequencies  between  1.2  and  14 
GHz,  a  restriction  which  is  not  important  for  our  present 
purposes. 

The  limited  field  of  view  afforded  by  the  use  of  large 
antennas  (27  and  40  m  diameter  paraboloids)  made  it  possible 
(at  all  but  the  lowest  frequencies)  to  restrict  the  field  of 
view  to  a  single  active  region.  The  observing  procedure 
consisted  of  sequencing  among  one  or  more  active  regions 
with  dwell  times  of  10  to  90  minutes  on  each.  Several  times 
during  the  day,  a  calibrator  source  (3C84  or  3C273)  was 
observed  tor  about  20  minutes. 

The  present  report  will  be  confined  to  the  first  3  of  6 
consecutive  days  of  observations  of  AR4549.  The  H-alpha 
appearance  of  this  region  on  the  first  day  (July  26,  1984) 
is  shown  in  Figure  2  where  it  is  seen  to  be  dominated  by  a 
single  sunspot.  Between  July  26  and  28,  the  spot 
progressively  weakened.  We  shall  be  interested  therefore 
not  only  in  the  coronal  field  distribution,  but  also  in  how 
this  distribution  evolved  as  the  spot  disappeared. 


III.  DATA  ANALYSIS 


The  fundamental  calibration  required  is  that  of  the 
interferometer  amplitudes  as  a  function  of  polarization, 
baseline  and  frequency.  This  calibration  was  based  on 
extended  observations  of  3C373  and  3C84.  Solar  amplitudes 
were  then  calibrated  by  direct  comparison  with  the 
amplitudes  observed  on  these  calibrators.  Noise  sources 
internal  to  each  receiver  were  used  to  correct  for  small  day 
to  day  variations  in  gain.  The  resulting  calibrations  are 
estimated  to  be  good  to  about  10%,  with  a  principal  source 
of  uncertainty  being  our  imperfect  knowledge  of  the 
calibrator  spectra. 

During  the  course  of  each  day,  the  observed  amplitudes 
and  phases  varied  as  the  rotation  of  the  earth  changed  the 
projected  baseline  (viz,  the  length  and  orientation  of  the 
baseline  as  seen  from  the  perspective  of  the  source).  These 
variations,  which  had  the  expected  qualitative  form,  can  be 
exploited  to  make  a  crude  map  of  the  source  at  each 
frequency.  Such  mapping,  however,  is  beyond  the  scope  of 
this  report.  Therefore  we  further  restrict  the  analysis  to 
data  acquired  in  a  "20  minute  period  each  morning,  a 
restriction  that  will  still  enable  us  achieve  our  objective 
of  measuring  the  size  and  brightness  temperature  spectra. 

The  use  or  a  common  time  of  day  minimizes  potential 
systematic  errors  introduced  by  the  neglect  of  the  details 
of  the  source  structure. 

The  fundamental  data  reduction  task  consisted  of 
automated  editing  of  the  data  to  remove  periods  of 
phase- lock  anomalies  or  marginal  antenna  tracking  and 
intervals  during  which  the  receivers  were  changing  frequency 
or  polarization.  The  data  were  then  sorted  by 
frequency/polarization,  and  corrections  made  for  the 
calculated  effects  of  baseline  errors,  earth  and  solar 
rotation.  Amplitudes  and  phases  were  then  vector-averaged 
over  the  chosen  ~20  minute  period,  and  the  amplitudes 
calibrated  as  discussed  above. 

This  procedure  yielded  fully  calibrated  amplitudes  and 
partially  calibrated  phases  on  each  of  3  days  for  each  of  3 
baselines  at  56  frequencies  in  both  right  and  left  circular 
polarization.  The  optical  data  showed  a  region  dominated  by 
a  single  sunspot.  Therefore  a  single  microwave 
gyroresonance  source  was  expected  as  well.  This  hypothesis 
was  tested  by  forming  the  closure  phase  (Pearson  and 
Readhead,  1984)  at  each  frequency  and  polarization.  Closure 
phase  can  be  defined  as  the  sum  of  the  three  phases  observed 
on  each  baseline.  Such  a  closure  phase  has  the  property 
that  its  value  is  independent  of  antenna-based  calibration 
errors  (viz,  calibrated  phases  need  not  be  used  to  form  the 
sum) .  Its  usefulness  here  is  derived  from  the  fact  that  the 
closure  phase  should  be  zero  for  the  expected  small,  single, 
symmetrical  source.  For  the  more  complex  sources  typical  of 


circumstances  where  free-free  opacity  becomes  important, 
this  would  not  necessarily  be  the  case.  Thus  at  each 
frequency,  the  closure  phase  represented  a  convenient, 
objective  parameter  which  provides  a  measure  of  the 
importance  of  gyroresonance  opacity.  Frequencies  and 
polarizations  for  which  the  closure  phase  did  not  come 
within  +-  l  radian  of  the  expected  value  were  excluded  from 
further  analysis. 

To  determine  the  size  of  the  source  at  each 
frequency/polarization,  we  assume  a  Gaussian  profile  and 
note  the  dependence  of  amplitude,  A  [Jy]  ,  on  projected 
baseline,  L  [km],  for  a  Gaussian  source  of  FWHM  diameter, 

D  [arcsec] ,  observed  at  frequency,  N)  1GHz]  : 
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where  A  [Jy]  is  the  limiting  amplitude  at  L  -  0, 
corresponding  to  an  unresolved  source.  Thus  a  plot  of  In (A) 
vs  L*  should  be  a  straight  line  whose  slope  yields  the  value 
of  D,  the  source  diameter.  A  typical  plot  of  this  type, 
shown  in  Figure  3,  shows  that  the  diameter  can  be  measured 
to  an  accuracy  of  about  5%.  It  might  be  noted  that  a 
principal  source  of  calibration  uncertainty,  the  flux  of  the 
calibrator,  does  not  affect  the  measurement  of  this  slope 
since  all  three  points  would  be  affected  by  the  same  error 
factor. 


Having  measured  the  diameter  and  AQ,  the  corresponding 
brightness  temperature  is  given  by  the  ratio  of  source  flux 
to  area,  using  the  relation. 


Tb  =  3.5*10*  AQ  (DW 


-2 


(3) 


The  resulting  size  and  brightness  temperature  spectra 
for  July  26  are  illustrated  in  Figure  4.  Figure  5  shows  the 
evolution  of  the  these  quantities  over  three  successive 
days. 


IV.  RESULTS 

Considering  first  the  overall  features  of  the  spectra 
shown  in  Figure  4,  for  the  LCP  brightness  temperature  we 
note  that  coronal  values  above  10*  K  are  maintained  up  to 
~1 2  GHz,  where  the  temperatures  drop  rapidly  to  a  few  times 
10*  K.  This  is  consistent  with  our  expectation  that  the 
brightness  temperature  of  the  source  would  drop  rapidly  at 
the  "transition"  frequency  where  the  corresponding  isogauss 
shell  was  tangent  to  the  transition  zone.  We  also  note  that 
at  lower  frequencies,  the  source  size  decreases  as  frequency 
increases.  The  relatively  constant  value  of  source  size 
above  "10  GHz  is  not  understood  at  present,  but  could  be 
related  to  spatial  fragmentation  of  high  frequency  sources 


such  as  ooserved  by  Lang  et  al.(1983)  at  15  GHz. 

Comparing  the  RCP  and  LCP  brightness  temperature 
spectra,  we  note  that  the  same  qualitative  behavior  is 
observed  in  RCP  as  in  LCP.  Of  particular  interest  is  the 
fact  that  shifting  the  RCP  spectrum  up  in  frequency  by  a 
factor  of  3/2  provides  a  reasonable  agreement  with  the  LCP 
brightness  temperature  spectrum.  This  is  consistent  with 
the  cnoice  or  the  2nd  and  3rd  harmonics  as  the  relevant 
values  for  the  ordinary  and  extraordinary  modes 
respectively. 

A  scaled  comparison  of  the  RCP  and  LCP  size  spectra 
might  also  be  expected  to  behave  similarly.  Although  they 
are  in  fact  consistent  in  terms  of  the  size  and  scaled 
frequency  at  which  the  size  levels  off,  the  RCP  results  at 
lower  frequencies  seem  significantly  higher  than  expected  on 
the  basis  of  the  LCP  spectra. 

An  unexpected  feature  noted  in  the  brightness 
temperature  spectra  below  the  transition  frequency  is  the 
apparent  systematic  increase  from  2  to  4*10*  K  as  frequency 
increases.  The  significance  of  this  feature  in  terms  of  the 
physics  of  the  coronal  fields  is  not  clear. 

The  evolution  of  the  coronal  field  as  the  sunspot 
disappeared  is  illustrated  in  Figure  5.  Assuming  that  the 
3rd  harmonic  is  the  relevant  one  in  the  extraordinary  mode, 
from  the  brightness  temperature  spectra,  it  is  apparent  that 
the  maximum  field  at  coronal  levels  decreased  from  about 
1300  to  1000  to  600  g  over  a  two  day  period.  The  relatively 
modest  change  in  the  size  spectra  at  the  lower  frequencies 
over  this  period  is  also  significant.  It  may  suggest  that 
the  weaker  fields  were  not  affected  by  the  spot's 
disappearance  to  the  same  extent  as  were  the  stronger 
f ields. 

Another  interesting  facet  to  the  field  evolution  was 
the  systematic  decline  in  brightness  temperature  associated 
with  the  strongest  fields  on  each  day,  from  4  to  2  to 
1*10*  K. 


V.  DISCUSSION 

The  observations  presented  above  represent  preliminary 
analyses  or  the  first  interferometric  size  and  brightness 
temperature  spectra  of  an  active  region.  The  region 
selected  for  this  analysis  was  deliberately  chosen  to  be  a 
simple  one,  a  single  isolated  sunspot.  In  broad  terms,  we 
have  seen  that  the  spectra  matched  expectations  based  on 
straightforward  application  of  gyroresonance  effects,  as 
exemplified  in  terms  of  magnetic  shells.  We  have  also  seen 
that  there  are  some  aspects  of  the  spectra  that  are  not  yet 
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understood.  With  more  sophisticated  analyses  and 
interpretations,  they  too  may  yield  their  secrets. 

At  this  preliminary  stage  of  the  analysis,  it  is 
inappropriate  to  consider  detailed  interpretations  of  the 
data  in  terms  of  the  solar  coronal  fields.  Nevertheless,  it 
may  be  worthwhile  to  review  the  strengths  and  weaknesses  of 
microwave  spectroscopy  (as  we  now  understand  them)  for  the 
measurement  of  such  fields. 

Unlike  photospheric  magnetograms  which  measure  magnetic 
FLUX,  (averaged  over  the  seeing-dependent  instrumental 
resolution),  microwave  spectroscopy  measures  magnetic  FIELD 
in  the  corona.  (Thus  the  present  data  are  inconsistent  with 
a  model  in  which  much  of  the  coronal  field  is  concentrated 
in  small  bundles  of  intense  field  as  has  been  suggested  at 
photospheric  levels.)  Through  equation  1,  the  field 
measurements  are  inherently  calibrated  so  that  long  term 
comparisons  are  quite  feasible.  Measurements  such  as  these 
can  be  made  in  a  few  seconds,  so  that  evolution  on  short 
time  scales  can  be  studied  also.  Full  calibration  of  the 
phases  will  enable  the  location  of  the  fields  to  be 
determined  interf erometrically  to  arcsecond  accuracy. 

The  use  of  microwave  spectroscopy  for  the  measurement 
of  coronal  fields  has  some  clear  weaknesses.  As  our  model 
suggests,  the  magnetic  measurements  represent  emission 
coming  from  a  range  of  heights,  characterized  by 
temperature,  not  distance  above  the  photosphere. 

Furthermore,  there  are  potential  sources  of  ambiguity 
related  to  the  specification  of  the  appropriate  harmonic, 
particularly  for  magnetic  fields  directed  close  to  the  line 
of  sight.  Perhaps  the  most  serious  weakness,  however,  is 
that  with  a  three-element  interferometer,  there  are  serious 
limitations  as  to  the  morphological  complexity  of  active 
regions  which  can  be  effectively  studied.  This  problem 
cannot  be  resolved  without  the  addition  of  more  antennas. 

In  the  near  future  with  presently  available 
instrumentation,  we  can  anticipate  the  achievement  of 
one-dimensional  multifrequency  synthesis  maps  of  simple 
unipolar  or  bipolar  regions;  more  sophisticated  solar 
interpretations  and  modelling  comparisons;  studies  of  active 
region  growth  as  well  as  decay;  measurements  of  height 
dependence  of  magnetic  fields  (from  limb  observations)  and 
application  to  studies  of  the  magnetic  fields  of  flaring 
loops. 

In  summary,  the  present  analysis  suggests  the  potential 
of  microwave  spectroscopy  to  provide  quantitative 
measurements  of  the  magnetic  field  distribution  in  the 
active  solai  corona. 
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FIGURE  CAPTIONS 


Figure  Is  Schematic  illustration  of  magnetic  field 
lines  above  an  isolated  sunspot.  The  isogauss  contours  for 
450 ,  600  and  900  gauss  are  shown.  The  4th,  3rd  and  2nd 
harmonics  of  the  local  gyrof requency  at  these  values 
respectively  are  equal  to  5  GHz. 

Figure  2:  H-alpha  filtergram  from  Big  Bear  Solar 
Observatory  for  AR4549  on  26  July  1984  when  the  region  was 
located  at  S17E07. 

Figure  3:  A  typical  plot  showing  the  dependence  of 
observed  amplitude  (circles)  on  projected  baseline. 

Gaussian  sources  of  different  FWHM  diameters  would  yield 
lines  with  the  slopes  as  indicated.  In  this  case  the  source 
FWHM  diameter  is  8.8  arcseconds. 

Figure  4:  Size  (top)  and  brightness  temperature 
(bottom)  spectra  in  LCP  and  RCP  respectively,  as  observed  on 
26  July  1984. 

Figure  5:  Evolution  of  the  LCP  (extraordinary  mode) 
size  (top)  and  brightness  temperature  (bottom)  spectra  as 
the  sunspot  disappeared.  Data  acquired  on  July  26,  27  and 
28  are  represented  by  0,  1  and  2  respectively.  With  the 
assumption  that  the  3rd  is  the  highest  optically  thick 
harmonic,  the  frequency  scale  can  be  interpreted  in  terms  of 
the  magnetic  field  scale  inserted  at  the  top. 
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